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Fig. 3. Fraction F of power radiated in backward direction for various 
antenna configuration and ground plane. @, is angle of maximum 
radiation from ground plane. 


The discrepaney between the measured and computed patterns 
in the backward direction may be attributed to two factors. Firstly, 
the diffraction field from the orthogonal edges have been neglected 
in computation but in practice it is large enough to influence the 
small field in the backward direction. Secondly, there is a small 
error in centering the dipole and the ground plane is not perfectly 
plane. One interesting point is that the shadowing by the ground 
plane is relatively insensitive to its width when it exceeds 5x. 
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Impedance Comparisons for the Asymmetrically Driven 
Thin Cylindrical Antenna 


JOSEPH H. PIERLUISSI 


Abstract—The results from various approximate methods of 
solution to three basic equations for the impedance of the asym- 
metrically driven perfectly conducting thin cylindrical antenna are 
briefly discussed and compared. Generally good agreement is found 
between method-of-moment solutions and the three-term King-Wu 
theory. Largest disagreement with those approximate methods is 
found with the Shen-Wu-King solution and the joining procedure. 


INTRODUCTION 


Although the radiation characteristics of the asymmetrically 
driven thin cylindrical antenna have been known from experiments 
for several years, rigorous theoretical studies of the problem were 
most probably instigated by Hallén’s analysis [1] of the center- 
driven case. For a perfectly conducting thin eylinder of radius a 
and length Z extending from z = —h, to 2 = h, and excited by a 
Dirac delta voltage source of magnitude Vo at z = 0, the basic 
equations adopted by succeeding workers to solve for the current 
T over the antenna surface may be written as: 


Vo +70Az = Voi(z2) (1) 
3 f# 
—+#)A, = Vos 2 
(e+ ) 03 (z) (2) 
PruVo 


A, = Cooskz +D sin kz — sink|z|. (3) 


No 


Here, the potentials A, and ¢ are the solutions to the Helmholtz 
equations, & is the propagation constant at a frequency w, «9 and 
mo are the medium characteristic parameters, and C and D are con- 
stants to be determined from the end-condition I (ti) = I(—h.) = 0. 
Clearly, (3) follows from (2) by integration, and (2) follows from 
(1) with the use of the Lorentz gauge. For solution, the approximate 
kernel is used in connection with A, and ¢. 


NuMERICAL SOLUTIONS 


The first practical solution to (3) was obtained by King and 
Wu [2] using an iterative procedure which lead to a three-term 
model for the current. Impedance calculations using the model 
were made by Harrison et al. [3], samples of which are shown under 
column (1) in Table I, The values missing on this table were not 
included in a similar table given in [8], in spite of the fact that the 
model is valid for the resonant lengths Z = 0.5 and Z = 1A. A 
much older and less useful iterative solution by King [4] resulted 
in an “averaging” method for obtaining the impedance of an asym- 
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cs 


CONDUCTANCE G (MILLIMHOS) 


TABLE I 


ImpepANCE VALUES as PROVIDED BY Various MODELS 


Total Arm 
Length Lengths 
L hy hy 

0.29 0.18 0.11 
9.29 0.21 0.08 
9.35 0.25 0.10 
0.35 0.29 0.06 
9.35 0.32 0.03 
9.48 0.26 0.22 
0.48 0.27 0.21 
0,48 0.32 0.16 
0.48 0.34 0.14 
9.48 0.37 0.11 
9.48 0.38 0.10 
3.48 0.40 0.08 
9.48 0.42 0.06 
9.48 0.45 0.03 
50 0,25 0.25 
9,50 0.30 0.20 
0.50 0.32 0,18 
0.5 0.35 9.15 
9.50 0.37 0.13 
0.50 0.40 0.10 
0.50 0.42 0.08 
9.50 0.45 0.05 
0.50 0.48 0.02 
0.80 0.48 0.32 
0.80 0.57 0.23 
0.80 0.64 0.16 
9.80 0.69 0.11 
9,80 9.72 0,08 
0.96 0.50 0.46 
0.96 0.51 0.45 
0.96 0.54 0,42 
9.96 0.58 0.38 
0.96 0.61 0.35 
0.96 0.64 0.32 
0.58 0.57 0.41 
0,98 0.73 0.25 
1.00 0.50 0.50 
1.00 0.55 0.45 
1.00 0.60 0.40 
1.06 0.70 0.30 
1.1) 0.57 0.54 
1.11. 0,59 0.52 
1.11 0,60 0.51 
1.11 0,62 0.49 
4.11 9.63 0.48 


—— SINUSOIDAL SUM 
--- JOINING PROCEDURE 
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(1) 
King-Wu | 


18,74-j 330.13 
19.26-4397.35 
34 .30-4256.50 
39.06-j365.41 
70.99+56.82 
73.43+j6.42 
92.59+53.23 
105.05+j0.54 
147.72-j12.69 
184, 23-j28.00 
238.05-4j61.97 
316 .09-4132.47 


680. 35-5290.17 
75.07-j 240.26 
51,19-4 209.98 
50.00-4 246.29 
42.19-}320.74 

801,56~3337.07 

619 .96-}219.40 

362.93-474.37 

191, 28-4 31.68 

132.40-59.13 
98.50-316.58 

287 .59+57.18 
89, 64+422.3% 


159 .45-j404.82 
196 .54-4425.18 
268,.11-4441.54 
370.05-3416.26 
440.57-3293.54 


ORIVING POINT AT 3/16 L FROM CENTER 


Input Impedance of Asymetrically-Driven Antennas 
(2) (3) (4) 


Shen- 
Wu-King 


25.77-3341.80 
24,.41-5407.35 
3842-5249 .66 
40.09-4346.01 
24.03-4598.98 
68.2444 21.06 
70. 324421. 31 
87.534422.79 
98.86+}23.20 
136.86+321.52 
170.33+519.84 
219 .44-43.54 
295.50~-429.53 
568.10-4 327.28 
77.184453.95 
86.124+)57.71 
99. 294+462.66 
122,51+469.83 
163. 734378.33 
240,73+482.87 
391.714+447.67 
651.06-4170.18 
746 .94-5862.44 
743.45-5227.28 
113.91-5259.76 
60.84-5213.67 
48.27-}230.51 
42,51-4275.12 
812.17-j288.62 
682.92-4174.36 
372,64-525.53 
217.09+42. 37 
145.40+}6.10 
108.86+45.58 
325.07+449 .80 
86.69+427.99 


669, 10-4538.16 
6A5.01-485.91 


£278 ,.894+5103.02 


106 .49+461.86 

268.16-5454.05 
307 .58-4459.10 
380 .87-j456.44 
490,18-3421.52 
612.03-j316.62 


2 LN(L/a)= 10 


Popovic- 
Surutka 


17.6-4318.8 
18.2-4377.4 
334-4244 .9 
37 .3-5336.2 
77544701 
79.9+46.8 
99.24+53.7 
111.7+j0.8 
153.3-413.2 
187,3-430.0 
234.9-}62.6 
298.0-j217.4 


492.0-4369.0 
89.8-4235.6 
51,3-j198.2 
42.9-5213.1 
41,8-j251.7 


Polynomial 
(3-terms) 


16.90-4301.46 
17.50-5356.56 
32,08-j229.87 
35.91-4313.75 
39.42-5527.94 
74 .02449.18 
76.294+49 .02 
94.974+56.86 
107.19+54,66 
147.47-46.96 
181.25-421.90 
228.93~j52.10 
292.92-4114.55 
391. 36-4469 30 
86. 724540.05 
96.924441.96 
111.96+443.99 
138.364)45.41 
184 .494441.76 
265.69+416.36 
392.44-488.58 
466 .08-4377,10 
299.43-5737.71 
442.81-5368.53 
79 .00-j217.11 
45.13-4182.04 
36.83-4195,56 
34.19-4229.57 


(5) 
Sinusoidal 
(3-terms) 


16 .92-4303,58 
17.54~4358.99 
31.99- 4231.54 
3584-5 315.77 
39 .59-5531.00 
72.00457.99 
75,164+47.81 
93.114]5.54 

104,774+43.35 

143.05-47.75 
175.06-4321.66 
220.45-449.29 
282.61-3105.95 

400. 75-4442.07 
85.314+438.78 
95.07+]40.63 

109.414+442.62 

134.464444.19 

178.014541.52 

254.88+520.40 

380.17-j70.07 
480.81-5342.41 
340.25~3737,11 

490.72-5 368.94 
83.23-4231.54 
45.49-4192.96 
36.01-4205.89 
32.72-4 240.22 


470.2-3462.2 416.49-3435.55 454,42-j441.46 
447,8-3345.2 414.30-4347.34 437.56-5347.95 
320.6-4140.0 317.82-j135.98 316.01-j141.94 


208.4-358.0 
144.7-529.5 
110.0-418.1 
332.9-444.1 

93.44428.5 


210.48-349.79 
147.75-j21.99 
113.19-J12.01 
309.77-§59.05 
96.744414.79 
287.21-4466.13 
422.B2~-4282.43 
293.76+48.69 
122.014+442.23 


206. 38-559.99 
144.40-431.98 
110, 20-j21,12 
302.15~j60,94 
91.68446.90 


322.13-4488.50 
442.80~J268.23 
283.28+4 7.86 
116.674434.81 


161.7-4364.5 134.82-j339.49 149.04-j362.42 
183,9-4372.5 153.38-4347.30 170.28-3369.95 
225.8-4382.2 188.24-4357.72 210.00-5378.94 
291.8-5385.0 243.85~j364.06 272.54~j380.83 
385.7-j361.9 322,63-j352,39 358.14-4358.77 


-—6 
wm 
fe] 
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—— SINUSOIDAL SUM 


--- JOINING PROCEDURE 
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10 
(6) 


Joining 


(7) 


Joining 


(Sinusoidal) (King-Middle ton) 


18.53-j297.37 
24.45~4 329.54 
47.49-5196.90 
71.61-j229-46 
115.01-j300.83 
74. 18+}8.81 
80.024511,09 
128.984) 25.27 
161, 40+] 30.27 
266.074525.11 
346.45-51.03 
442,44-563.66 
534.36-4181.55 
561.89- 5539.70 
85.31+538.78 
109 .90+347.63 
146.414} 56.76 
209.794+462.47 
312.77+}46.06 
454. 76-339. 32 
565.83-4248.60 
529.19-j515.52 
387. 78-678. 24 
393.43-4173.65 
97.63-5 202.81 
47.25-j196,18 
38.10-4195.00 
42,59-5196.98 
466.48-J490.74 
481, 76-]453.44 
493. 20-}281,07 
396.74~4104,71 
264.71-326.23 
168. 35~-j10.03 
458.88-4227, 31 
99,.83+46.62 


322.13-3488,50 
393.57-j415.89 
423,88-j150.39 
135.23+j34.21 

145,82-j 359.51 
156.41-3358.85 
174.87-j358.68 
201.78-j351.12 
237.54-j339.03 


20.11~4 306.89 
25.95-4341,10 
48.595 204.50 
7158-5245, 72 
74.824+48.87 
80.17+j11.11 

127,12+425.23 

159.074+430.11 

267.6944 22.39 

351. 74-39.42 

450,99-484.91 

538,89-4224,95 
352.83-4159,19 
97,65~4186.95 
53,59-4199.77 
43. 46-4198,05 
46.47-5210,84 
418.92-4472.42 
447.24-4430.4] 
481.40-j251,83 
382.59-569,.08 

247 .32-33.65 

158,58-jl.11 

&48,87~]221.37 
96. 30+424.63 


137.19-4320.67 
145.78-4321.46 
161.18-5321.98 
184.90-5320.60 
218.71-4314.22 


ORIVING-POINT AT 3/16L FROM CENTER 


2LN(L/Za)®* 10 


(b) 
(a) Conductance of asymmetrically driven dipole. (b) Susceptance of asymmetrically driven dipole. 
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metrically driven antenna from the respective impedances of two 
eenter-driven antennas, having arm-lengths equal to each one of the 
arm-lengths of the asymmetric case. The method was called the 
‘Joining procedure” in [8] and used, together with the well-known 
King—Middleton model [5] for center-driven antennas, to generate 
the values listed in Table I under column (7). 

If the first two terms to the right of (3) are discarded and the 
proper integration limits placed on A., it results in an integral 
equation for the current on an infinitely long antenna for which an 
exact integral solution is available [6]. The integral was evaluated 
in approximate analytical form by Shen et al. [7], and the results 
used to propose a three-term traveling-wave expression for the 
current. The model was used by the present author to calculate 
the impedance values given in Table I under column (2), as well 
as to generate the admittance curves of Fig. 1. 

Popovic [8] solved (8) with the method of moments assuming 
as basis function a current distribution in the form of polynomials 
in 2, continuous over each antenna arm and forced to be equal at 
the driving-point. Using third-order polynomials and point-matching 
the present author computed the impedance values shown in Table 
I and under column (4). The center-driven ease of (3) was studied 
by Neff eé ai. [9] using a three-term sinusoidal series for the current, 
defined over the entire antenna length, together with the method 
of moments and point-matching. The present author adopted their 
sinusoidal series for the current J (z) and defined it over each antenna 
arm, that is, 


{» 
SD A,sin — (i —2), O<2<mt 
n=l 2h 
I(z) = ; 
D Ba sin 5 (ha +2), —-h<2<0 (4) 


to compute with the same procedure the impedance values shown 
in Table I under column (5) and the admittance curves of Fig. 1 
labeled “sinusoidal sum.” Furthermore, with the center-driven 
case of (3) and (4) and the joining procedure column (6) in Table 
I was prepared, as well as the admittance curves of Fig. 1 labeled 
‘§Soining procedure.” 

Equation (2) was first used by Harrington and Mautz [10] in 
connection with the method of moments and point-matching, 
assuming a triangular- (as well as pulse-) subsectional basis function 
for the current. Unfortunately, the impedance tables referenced in 
[10] do not include the case Q = 2In (L/a) = 10 used in the 
preparation of Table I and Fig. 1. Popovic and Surutka [11] 
generalized the two-term variational procedure of Storer [12] to 
analyze (2). Their results are listed in Table I under column (3). 
Although values for the resonant lengths are missing in this column, 
the method is good for any kh, or khe less than 37/2. A detail de- 
velopment of (1) for use with the method of moments and pulse 
basis-function expansion may be found in the book by Harrington 


[13]. 
Discussion 


A close examination of Table I reveals good agreement between 
the impedance values under columns (3), (4), and (5), all of which 
may be considered to be method-of-moments solutions. The agree- 
ment between these and the King-Wu values in column (1) is 
reasonable as a whole. The reactance components from the Shen- 
Wu-King model [column (2)], however, disagree significantly in 
some eases from the values given by the other two basie methods. 
The disagreement with a method-of-moments calculation (i.e., 
sinusoidal series) may be more clearly observed in Fig. 1. The 
poorer agreement between the previous results and the joining 
procedure, as seen both in column (6) and (7) and in Fig. 1, is to 
be expected since the latter is only a very approximate method. 
The curve labeled “sinusoidal sum” in Fig. 1 is in excellent agree- 
ment with [10, fig. 2]. Impedance calculations were made by the 
present author using (1) but they did not show significant difference 
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with the other method-of-moments values. For instance, for the 
half-wave dipole with L/2a = 200 and 31 subsections, equation (1) 
gave an impedance of 85.26 + 743.21, while Harrington and Mautz 
listed 86.4 + 745.4 in [10]. 

All of the models described can provide impedance values for 
the center-driven antenna, but none can treat the end-driven case. 
Except for the joining and variational methods all can yield the 
current distribution along the antenna. All the current-distribution 
models excepting the model of Shen et al. make the initial assumption 
that the current vanishes at the antenna ends. By far the one 
requiring the least computational time for a given antenna is the 
model of Shen et al. This is in view of the fact that to use the at- 
tractive joining procedure advantageously one requires a non- 
existent simple model for the center-driven case. It is recalled that 
a truneated-Fourier series fits the function described in a least- 
squares sense, and as such the accuracy of the representation 
depends on the location of the sampling points used in the analysis. 
For the reported computations with (4) the sampling points needed 
in the determination of A, and B, were evenly distributed along 
each antenna arm, with one sample taken always at the driving- 
point. A reasonable criterion to use with (4) is that the number of 
terms chosen be so that one of them have the period of the sinus- 
oidal-current assumption for the specified length. Hence, the 
longest antenna that can be treated with three terms of (4) is the 
center-driven ease of L = 1.50x. This value is also the limiting length 
for the King~Wu model. Based on their [7, fig. 3], Shen eé al. in- 
dieated that their model was good only for the cases where the 
driving point was farther than about 0.15A from either antenna end. 
Calculations made by this author consistently gave reasonable 
values (see Table I) whenever L > 0.3 for a wide range of driving- 
point positions, the precise length being dependent on the antenna 
radius. 

The impedance values in Table I are intended primarily to 
extend to seven cases the two-case comparison of Harrison e¢ al. 
[3], thus updating their results with new models and with the 
presently used method-of-moment methods. Perhaps the only 
intelligent observation that can be made from the results given 
here, is that a comprehensive theoretical and experimental study is 
greatly needed to establish the accuracy of the several methods 
available for the impedance (and current distribution) of the asym- 
metrically-driven linear antenna. 
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Synthesis of Broad-Band Arrays with Arbitrary 
Frequency-Independent Elements 


FRANC BRATKOVIG 


Abstract—The application of dynamic programming to the syn- 
thesis of broad-band antenna arrays with arbitrary frequency- 
independent elements is investigated. Four criteria are proposed 
for the optimization of the array. By employing these criteria it is 
shown that the optimum array is the one with monotonic spacings. 
For such arrays the principle of optimality is proved and the limits, 
within which dynamic programming is applied, are analytically 
calculated. A 35-element array with ratio of upper-to-lower fre- 
quencies equal to 3 is synthesized by this method. 


INTRODUCTION 


Grating lobes occur in the radiation pattern of a broad-band 
antenna array. Suitable spacing of elements reduces the grating 
lobes. The optimum configuration of element spacings in the broad- 
band sense is the one that produces minimal peak magnitude of 
the radiation pattern in the visible range outside the main lobe in 
the whole frequency band. 

Ishimaru [1] transformed the radiation pattern into Poisson’s 
sum. Ishimaru and Chen [2] reported that the radiated power of 
each grating lobe is equal to the radiated power of the main lobe, 
regardless of the spacing and the excitation of the elements. Applying 
stationary phase approximation Chow [8] synthesized the op- 
timum nonuniform array of isotropic radiators. He also investigated 
the effect of the element pattern, and in turn reduced the first 
grating lobe of a narrow-band thinned array with an element pattern 
having a minimum coinciding with the first grating lobe. 

A basically different approach to the synthesis of the optimum 
nonuniform arrays was purely numerical. Skoinik et al. [4] applied 
dynamic programming to nonuniform arrays, however, without prov- 
ing the applicability of the principle of optimality. Since their method 
involved a wide range of possible spaces, quantization of the latter 
was coarse, to maintain computer time within reasonable limits. 

In the present communication an analytical method is combined 
with dynamic programming in order to optimize broad-band arrays 
with arbitrary frequency-independent elements. Dynamic pro- 
gramming is applied to the array within analytically calculated 
boundary values. It is proved that the principle of optimality 
applies. 


ForMULATION OF PROBLEM 


A linear array (Fig. 1) of N nonuniformly spaced elements with 
radiation pattern F’.(@) can be described by two vectors, the com- 
plex current vector J = (11,Io,+++,fw), and the real position vector 
X = (%1,%9,+°+,4y). The radiation pattern of the array is 


P.ixu) = Feu) Feu), ue (*) ain Gi) 
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Geometrical configuration of unsymmetrical nonuniform linear 


Fig. 1, 4 
antenna array with N elements. 


where the array factor is 


N 


F(ix,u) = 3) Inexp (jrau). (2) 


n=l 


The frequency-limit ratio of an antenna is the ratio between the 
highest and the lowest frequency at which the antenna still has the 
prescribed properties. In an antenna array the lowest frequency f; 
is limited by the minimum space smin, which must be no less than 
half the wavelength \;. The highest frequency fe is determined by 
the highest permissible boundary of the visible range wan. Thus 
the frequency-limit ratio is 


R, =f = Senta (3) 
fi © 


The frequency-limit ratio of broad-band arrays depends mainly 
on the array factor. By introducing Poisson’s summation formula, 
the array factor is transformed into 


mo(u) 
F(Ljnu) = 3. Pa (Tu) (4) 
meom(u) 
NU2 
Fy(i,aju) = f I(v) exp { j[ux(v) — 2mrv]} de. 
1/2 


Discrete variable n is substituted for by continuous variable v. 
Vectors J and x are substituted for by continuous functions I (v) and 
x(v) fulfilling conditions I(n) = I, and 2(n) = tn, n = 1,2,-+-,N. 

The term Fy represents the main lobe, whereas the other terms 
F, represent the grating lobes. The terms m,(u) and m(u) are 
the limits of the significant terms. The set of the spaces 
Sn = Inui — Ln, M = 1,2,°-+,N — 1, can be also substituted for by 
continuous space function s(v), which fulfills conditions s(n) = s,, 
n = 1,2,+++,N —1. Using the stationary phase method [5] and 
the approximation z’(v) = s(v) one finds the interval in which F,, 
is significant: 

Sue. 6) 


Smin 


2mr 


Smax 


The magnitude of the lobe within the interval is 


eta) Ea s=< ou). ©) 


Nm? | s'(v) 


The variable » depends on u because of the stationary phase 
condition 


us(v) — 2mmx = 0. (7) 


